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ABSTRACT 

Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. 
This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. This process is critical for 
the galaxy evolutions. We have investigated several aspects of the process as the timing and duration of the dilution and 
its correlation with the induced star formation. We analysed major (1:1) gas-rich interactions and mergers, spanning 
a range of initial orbital characteristics. Star formation and metal enrichment from SNe are included in our model. 
Our results show that the strongest trend is between the star formation rate and the dilution of the metals in the 
nuclear region; i.e., the more intense the central burst of star formation, the more the gas is diluted. This trend comes 
from strong inflows of relatively metal-poor gas from the outer regions of both disks, which fuels the intense star 
formation and lowers the overall metallicity for a time. The strong inflows happen on timescales of about 10 8 years 
or less (i.e., on an internal dynamical time of the disk in the simulations), and the most intense star formation and 
lowest gas phase metallicities are seen generally after the first pericentre passage. As the star formation proceeds and 
the merger advances, the dilution reduces and enrichment becomes dominant - ultimately increasing the metallicity of 
the circumnuclear gas to a level higher than the initial metallicities of the merging galaxies. The "fly-bys" - pairs that 
interact but do not merge - also cause some dilution. We even see some dilution early in the merger or in the "fly-bys" 
and thus do not observe a strong trend between the nuclear metallicities and separation in our simulations until the 
merger is well advanced. We also analyse the O and Fe enrichment of the ISM, and show that the evolution of the a/Fe 
ratios, as well as the dilution of the central gas metallicity, can be used as a clock for "dating" the interaction. 
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1. Introduction mergers and allow us to constrain the relative time scales 

of the different epochs (i.e., date the merger /interaction). 

Understanding how galaxy interactions and mergers shape If the Progenitor disks have strong gas metallicity gra- 
the ensemble population of galaxies in the local universe dients as observed m many galaxies m the local um- 
is crucial for our understanding of galaxy evolution. It has verse flShields| H22fl i Dmersteml 11223), ™e would expect 
been recognised for several decades that mergers may play that interaction-induced gas inflows will drive a noticeable 
an important role in th e evolution of galaxies, especially amount of gas into the central regions from the outer disk 
for the early types (e.g., iToomre I fl977h . While there have and thus initially lower the circumnuclear gas phase metal- 
been many observations and theoretical studies of "first- llclt y- Outflows of gas ejected m tidal tails may contribute 
! order" effects in interactions and mergers (such as morpho- to modifying the metallicity profile of the final merger rem- 
logical and kinetical evolutiona and star formation), much nant - Subsequent star formation (and outflows) may in- 
Ctf 1 less work has been done on 2 nd order effects such as the crease (decrease) the gas phase circumnuclear metallicity. 
evolution of metallicity ^he exac ^ timing of the inflows versus the initiation of in- 
tense star formation and subsequent metal enrichment may 
Recently, we have undertaken the study of the evolution provide robust constraints on the underlying physical mech- 
of metallicity in major dry (i.e., gas poor) mergers, show- anisms that determine the rate and relative timing of gas 
ing that they can lead to metalli city gradients in agree - flows an d star formation within the models. Comparing the 
ment with those of local ellipticals (jDi Matteo et al.ll2QQ9h . evolution of the metallicity and its relationship to both in- 
However, in gas-poor mergers, the evolution of the local flows and star formation will allow us to constrain the dis- 
galaxy metallicity is driven only by the mixing of stars dur- sipation time scales, the star formation evolution, the time 
ing the coalescence process. Gas can alter the properties of delay between the onset of inflows, sufficient to dilute the 
the final remnant significantly, affecting both the final gas metallicity, etc. 

phase and stellar metal abundance and ratios. In partic- In the past few years, a number of studies have in- 

ular, gas phase metallicity can play an important role in vestigated the role played by interactions and environ- 

constraining our understanding of the complex interaction me nt in determining the gas p hase metallicities of galax- 

between gas flows and star formation in interactions and ies. Donzelli & Pastoriza (2000) show that merging galax- 
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ies have on average higher excitation in their optical emis- 
sion line spectra than interacting pairs, and attribute 
this difference to low er gas metallicity in the mergers. 
Ma rquez et alJ (|2002f ) studied a sample of more than one 
hundred spiral galaxies, ranging from isolated to interac- 
tion with strong morphological distortions, and show that 
the [Nil] /Ha ratios, used as a metallicity indicator, in- 
dicate a clear trend from the metallicity to morphologi- 
cal type, with earlier type spirals showing higher ratios, 
while they found no trend wit h the s tatus of the inter- 
action. Recently, iKewlev et al.l ([20061 ) have derived the 
luminosity-metallicity relation for a sample of local galaxy 
pairs and compared it with that of nearby field galaxies. 
They found that pairs with small projected separations (s 
< 20 kpc h _1 ) have systematically lower metallicities than 
either isolated galaxies or pairs with larger separations, for 
a given luminosity. They also found a correlation between 
gas metallicity and burst strengths - all galaxies in their 
interacting sample with strong central bursts having close 
companions and metallicities lower than the comparable 
field g alaxies or pairs with wider separations. IRupke et al.l 
([2008t h in a study of strong interactions with high star for- 
mation rates (ultra-luminous infrared galaxies), found that 
the metal abundance in these intense starbursts is a factor 
of two lower than that of galaxies of comparable luminos- 
ity and mass. These results have been generally explained 
by gas inflows induced by the interaction, which dilute the 
pre-existing nuclear gas to produce a lower metallicity than 
that observed for wider separated pairs or isolated systems. 
More recently, iPeeples et all ([2008L I2009D have studied, re- 
spectively, a sample of low- mass, high- metallicity and high- 
mass, low- metallicity outliers from the mass- metallicity re- 
lation of star-forming galaxies from selected from the Sloan 
Digital Sky Survey (SDSS).They showed that the low-mass, 
high-metallicity outliers are usually isolated galaxies, with 
no evident companion or strong interactions. On the other 
hand, the high-mass, low-metallicity outlierts typically con- 
sist of systems that have high star formation rates and ev- 
idence for disturbed morphologies. 
However, ICooper et al.l 



([20081 ) found a strong 
metallicity-density relation for star forming galaxies 
in the local universe, with the more metal-rich galaxies 
apparent ly favouring region s of higher galaxy overdensity 
(see also lEllison et al.ll2009t ). They conclude that the dis- 
crepancy found with the outer studies (including those just 
discussed) is due to the fact that the number of close pairs 
(s < 100 kpc h _1 ) in the SDSS sample constitutes only a 
tiny fraction of the whole sample (less than 1% of galaxies). 
In this case, close pairs with low central metallicity could 
not contribute significantly to the scatter found in the 
mass-metallicity relation. Of course it can be difficult 
to discern trends in the gas phase metallicity during a 
merger when studying large samples of galaxies in SDSS 
as other factors may dominate the overall mass-metallicity 
relationship and its scatter. A confirmation that the 
mass-metallicity relation is affected by interactions only 
for close pairs showing signs of str o ng di sturbances has 
been found by iMichel-Dansac et~aT1 ([20 08), who pointed 
out that, in such pairs, the gas metallicity depends on the 
mass ratio of the two interacting systems: less massive 
members are systematically enriched, while a galaxy in 
interaction with a comparable stellar mass companion 
shows a metallicity lower than that of a galaxy in isolation. 
So while mergers and interactions alone may not drive the 



overall scatter in the mass-metallicity relationship, it may 
be a contributing factor. All these studies suggest that the 
dilution and enrichment of the ISM in the central regions 
of disk galaxies strongly depend on the exact timing of the 
different processes at play: gas inflows, interaction-driven 
star formation, gas consumption, feedback and subsequent 
enrichment. 

While a number of numerical studies have investigated 
the response of the gaseous component of galaxies during 
tidal interactions an d the subsequent star formation such 
interactions induce ([Iono et al.l [20041 : Mihos & Hernquist 
| l994atlSDringelll2000l:ICox et a l. 2006, 120081 : iKapferer et al. 
l2QQ5t iDi Matteo et al.ll2007ll2008D . little attention has been 
given to the detailed evolution of the metal content dur- 
ing galaxy encounters. Usi ng a galaxy pair ca talogue from 
cosmological simulations, iPerez et "all (|2006t ) have shown 
that the O/H abundance ratio in the central regions of 
close galaxy pairs (s < 50 kpc h _1 ) shows a lower level 
of enrichment than the mean O/H abundance ratio of a 
control sample, thus confirming the role played by gas in- 
flow in d iluting the metal co ntent of the nuclear regions. 
Recently, IRupke et al.l (|2010t ) have analysed simulations of 
major galaxy mergers, studying the dilution of the gas 
metallicity in the nuclear regions due to gas inflow. They 
found a dilution of about 0.1 - 0.3 dex, happening shortly 
after the first pericentre passage between the two galax- 
ies. However, their models do not include either star for- 
mation prescriptions or metal enrichment, so that while it 
has been possible to give predictions about the strength of 
the dilution, nothing is still known about the role played by 
interaction-induced star formation in the metal dilution and 
then subsequent enrichement. Moreover, the exact timing 
of the dilution peak and its correlation with the increase in 
the amplitude of the star formation have not been studied 
yet in any detail. 

2. Models and initial conditions 

Here we study mergers and flybys involving two massive 
Sbc galaxies, having a mass ratio 1:1. The Sbc galax- 
ies (hereafter called gSb) are composed of a spherical 
dark matter halo and a spherical bulge, rep resented by 
Plummer spheres ([Binnev fc Tremaind 119871 ) with total 
masses M H = 1.7 x 10 11 M and M B = 11.5 x 10 9 M 
and core radii th = 12 kpc and tb = 1 kpc respectively. 
The stellar disk is represented by a M iyamoto-Nagai den- 
sity profile ([Binnev fc TremaineTll987t ) with mass M* = 



40.6 x 10 9 Mq, vertical and radial scale lengths given re- 



5 kpc. The galaxies ini- 



0.2 M*, redistributed in a 



spectively by h* =0.5 kpc and a 
tially contain a gas mass M gas = 
Miyamoto-Nagai disk with vertical and radial scale lengths 
given respectively by h gas = 0.2 kpc and a gas = 6 kpc. For 
each pair of interacting galaxies, we performed 24 simula- 
tions, varying the galaxies orbital initial conditions (initial 
orbital energy E and angular momentum L) and taking into 
account both direct and retrograde orbits. The initial or- 
bital parameters (initial distance between the two galaxies, 
initial relative velocity, specific angular momentum, orbital 
energy and s pin) for the differen t runs are fully described 
m Table 7 of lChilin garian et all ([2010D . and we refer the 
reader to this paper for their complete description. 

We chose a reference frame with its origin at the 
barycentre of the system and x-y plane corresponding to the 
orbital plane. For each interacting pair, we have kept the 
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disk of one of the two galaxies in the orbital plane (i\ =0°), 
and varied the inclination i 2 of the companion disk, con- 
sidering: i 2 = 0°, i 2 = 45°, i 2 = 75°, and i 2 = 90°. The 
orbital angular momentum can be parallel (direct orbit) or 
anti-parallel (retrograde orbit) to the z-axis of the reference 
frame. 

All the simulations (96 in total) were run u si ng th e 
Tree-SPH code described in Se melin fc Combesl ([20021 ). 
Each galaxy contains a total of Ntot — 120000 particles, 
equally redistributed among gas, stars and dark matter. A 
Plummer potential is used to soften gravity at small scales, 
with constant softening lengths of e = 280 pc for all par- 
ticles. The gas is modelled as isothermal, with a tempera- 
ture T gas = 10 4 K. The assumption of an isothermal ISM 
is justified by the short cooling times in the gaseous disk at 
the resolution of our models, such that fluctuations in the 
gas temperature, those due to SN explosions for example, 
would be quicky radiated away (see the discussion in, e.g., 
iMihos fc Hernauistl[T996h . 

We refer the reader to ISemelin fc Combesl (|2002f ) for a 
detail description of and the tests perform ed to validate the 
code. Other tests have been presented in iDi Matteo et all 
(|2008[ ). The equations of motion are integrated using a 
leapfrog algorithm with a fixed time step of At = 5 x 10 5 yr. 



2.1. Star formation and metal enrichment 

The numerical recipes used to implement star formation 
and metal enrichment have been fully described elsewhere 
(|Di Matteo et al.l l2007: C hilingarian et al.ll201Q[ ). thus here 
we only summarise those features that are most directly 
relevant to our analysis. Star formation efficiency has been 
parametrised as: 



gas 



C x 



Pgas 



1/2 



(1) 



where M gas is the gas mass of a particle, and p gas the 
local gas density. The constant C = 0.3 pc 3 / 2 Mq 1 ^ 2 Gyr -1 
has been chosen such that the isolated disk galaxies form 
stars at an average rate of ~ 1 Mq yr -1 . Once the star 
formation recipe is defined, we apply it to gas particles, 
using t he hybrid method described in IMihos fc Hernquistl 
(|l994bl ). 

The effect of star formation on the surrounding in- 
terstellar medium has been implemented as follows. For 
each star- forming hybrid particle, we evaluate the frac- 
tion of stars f o rmed with masses > 8 Mq, adopting a 
iMiller fc Scalol (|l979l ) initial mass function (IMF), and 
we assume that stars above this mass threshold instanta- 
neously become super novae. Each of the SN leave behind 
a remnant of 1.4 Mq and releases their remaining mass to 
the surrounding ISM. The mass released also enriches the 
surrounding gas with metals. This is done assuming a yield 
y = M ret /M* = 0.02, where M ret is the total mass of all 
reprocessed metals and M* the total mass in stars. For each 
gas particle, the return of mass and metals is applied to the 
i — th neighbour gas particle, using a weight Wi based on the 
smoothing kernel. The metallicity is initially distributed in 
a gradient of the form: 



z m (R) = zo x 10 



-0.07 R 



(2) 



with R the particle distance from the galaxy centre 
and z = 3 z Q (jKennicutt et al.ll2003l : iMagrini et all 120071 : 



iLemasle et al.l 120 08), evolves with time, as star formation 
proceeds. 

Supernovae explosions also inject energy into the sur- 
rounding ISM. This is taken into account assuming that a 
fraction, e^ n , of the energy, Esn, released by a SNe goes 
into kinetic energy giving a radial kick to the velocities of 
neighbouring gas particles. The value of e/^ n has been cho- 
sen such that the total amount of kinetic energy received 
by a gas particle due to the contribution from all the sur- 
rounding neighbours, is less than 1 km s _1 , thus preventing 
a rapid growth of the vertical thickness of the gaseous disk. 
A more detailed description of the impl ementation of feed- 
back re cipes used in our code is g iven in lSemelin fc Combesl 
ti2002D : IChTlingarian et all (120101) . 



3. Results 

3.1. Gas Flows: Where and When 

It is known that close passages of galaxies can destabilise 
galaxy disks, producing asymmetries, such as bars, which 
are efficient mechanisms for driving gas into the circumnu- 
clear region of a galaxy, produc ing a burst of star formation 
(|lono et al.ll2004tTComb es 2008). If, prior to the interaction, 
the gas is distributed in the disk with a negative metallic- 
ity gradient, with the outer regions being relatively more 
metal poor than the inner ones as observed in local galaxies 
([Shiel ds 1990: lDinersteinlll996[ ). then a portion of this low- 
met allicity gas will fall into the nuclear regions, diluting the 
pre-existing gas with gas of lower metallicity. This dilution 
will last until the star formation, which is enhanced by the 
gas inflow, releases reprocessed metals to enrich the gas. 

This cycle (dilution and enrichment) is shown in Figs. 
[T] and [2] for some of our galaxy merger simulations. For the 
analysis of the effect of mergers on the metallicity of galax- 
ies, we have defined the metallicity dilution (often referred 
to as simply dilution) as z/zi SO < 1, where z is the metal- 
licity of a galaxy in the pair during the merger or flyby and 
Zi SO is the corresponding value for the same galaxy evolved 
isolated. Both z and Zi SO are measured in apertures of 1-2 
kpc and thus represent the evolution of the metallicity for 
the nuclear or circumnuclear regions of galaxies (which also 
corresponds reasonably well with the aperture sizes used to 
estimate the metallicities of galaxies; see discussion in the 

§i). 

We find that as soon as the two galaxies have undergone 
their first pericentre passage, an intense inflow of gas, last- 
ing for 300-400 Myr, takes place. A detailed discussion of 
the mechanism genera ting the gas inflows we o bserve in our 
simulations is given in IDi Matteo et al.l (|20Q7h . They show 
that each close passage of the two interacting galaxies is 
accompanied by the amplification of stellar asymmetries - 
stellar bars. These stellar bars are responsible of removing 
angular momentum from gas, transferring it to the stellar 
component, and in regulating the subsequent inflow of gas 
in the inner kpc. We note that this result has already been 
discussed in several previous studies using simulations (see 
for example Mihos & Hernquist 1994, 1996). This mecha- 
nism does not require any stellar feedback to work, that 
is, the gas inflows are a natural consequence of the dy- 
namics of the interaction and the redistribution of the an- 
gular momentum that is induced dur ing the prog ression 
of the merger (|Di Matteo et all 120081 : iRuuke et all 1201 Ol ). 
This intense inflows can deposit infewl0 8 ys into the cen- 
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Fig. 1. The dynamics of the gas during the interaction and subsequent merger of two gSb galaxies for several orbits 
(shown with different colours). (First panel): Relative distance of the two galaxy centres versus time for some simulated 
orbits. (Second panel): Gas mass inside 2 kpc from one of the two galaxy centres versus time. (Third panel): Evolution of 
the star formation rate (SFR) relative to the rate for the corresponding initial galaxies evolved in isolation. (Fourth panel): 
Evolution of the central gas metallicity. The metallicity has been normalised to that of the initial galaxies used in the 
merger simulations evolved isolated. In the second, third, and fourth panel, squares and asterisks respresent respectively 
the time of the first pericentre passage and the time of final coalescence for the different orbits. 
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Fig. 2. Same as FigJTJ but for the simulations of "flybys". Note that simulations resulting in starbursts with higher 
intensities also have gas phase metallicities with higher dilution. Even in flybys, the star formation rate and dilution of 
the metallicity are related. 



tral (r < 2 kpc) regions ~ 2 — 3 x 10 9 Mq of gas, origi- 
nally located before the interaction at an average distance 
of 6 — 7 kpc from the galaxy centre. The amount and origin 
of the inflowing gas obviously depends on characteristics of 
the merger like orbital parameters of t he interaction and 
gas fraction (see iDi Matteo et al.l 120071 for a discussion). 
In all cases, this inflow of gas fuels the first increase in the 
star formation rate during the merging process. A second, 
sometimes even stronger peak of SF takes place in the final 
phases of the interaction when the two galaxies are close to 
coalescence. 

The time of strong enhancement in the SFR also coin- 
cides with the time of high-dilution. Typically, we find that 



the dilution lasts for 2 x 10 9 years, with half of the sample 
substaining a z/zi SO < 0.8 for less than 5 x 10 8 (Fig]3]). In 
the case of fly-bys, this dilution is still visible far after the 
pericentre passage, when the separation between the two 
galaxies is several hundreds of kpc. We emphasise that the 
gas redistribution and subsequent dilution takes place on 
scales of several kiloparsecs and is not due to the effects of 
stellar feedback which typically act over much smaller phys- 
ical scales. Moreover, significant star formation is initiated 
only after a substantial gas inflow has occurred. So feed- 
back from SNe explosion has little effect on the inflowing 
gas and certainly is unable to prevent substantial inflows of 
gas initially. 
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A possible concurrent physical mechanism for redis- 
tributing metals in the interstellar medium of the galaxy 
could be through gas diffusion, which is not included in our 
cod^E At scales of hundreds of parsecs (which corresponds 
roughly to our numerical resolution), the diffusion law can 
take the form of a heat equation, with a diffusion coeffi- 
cient D = 9.25 x 10 26 cm 2 s~ 1 (see iMart Inez- Serrano et aD 
120081) . This means that the typical timescale for diffusing 
metals from the outer regions of the disk (R =7kpc0) is 
Tdiff — R 2 /D > 10 10 yr. This timescale is significantly 
longer than the average inflow timescale which is of-order 
300-400 Myr. Thus, we can safely assert that the process of 
dilution of the central gas metallicity studied here is driven 
by interaction- driven inflows rather than metal diffusion. 

Once the gas inflow and the corresponding burst of star 
formation is over, the re-processed gas released by SNe ex- 
plosions starts to enrich the central metallicity, so that the 
z/zi SO ratio becomes greater than one (Figs. [T] and [2j fourth 
panel). We checked that during the interaction, the con- 
trast in the gas density between the nuclear region and the 
outer disk are sufficient to result in gas depletion times that 
are significantly lower than the times typically required to 
re-enrich the ISM in metals. This substantiates our inter- 
pretation that the dilution in the metallicity of gas in the 
nucleus is due to strong merger-driven gas inflows. 

3.2. Relationship Between the star formation and Metallicity 

We have seen that the dilution of the gas-phase metallicity 
takes place at approximately the same time as the strongest 
enhancement in the star formation. In general, we find that 
the peak of SF and the maximum dilution are roughly co- 
eval, with a delay between the two of 10 8 yr at most. This 
correlation between the starburst phase and the strong dilu- 
tion is clearly seen in the results of the simulations shown in 
FigHJ The maximum dilution is reached in the simulations 
at approximately the peak in the star formation rate. In 
this phase, the higher SFR, which is a sign of the strongest 
gas inflow are indeed associated with the strongest dilution 
as would be expected. We note that a strong correlation 
between metallicity and centra l burst strength has been ob- 
served bv lKewlev et al.l ([20061 ) in a sample of galaxy pairs. 
Specifically, they pointed out that all five galaxies in the 
pairs sample with strong central bursts have close com- 
panions and metallicities lower than the comparable field 
galaxies. 

Before the peak in the star formation rate, galaxies tend 
to populate a region of low to moderate SFRs and dilutions 
(typically 0.7 < z/zi SO < 1; see FigHJ). During the burst 
phase and peak of star formation rate, the eject a of the 
SNe starts to enrich the gas, and indeed, once the burst 
is over, galaxies have low star formation rates and high- 
metallicities (FigHJ). The first two phases of the relation 
between star formation and dilution are associated with 
before, during, and after the first pericentre passage and 
then close to final coalescence of the merging galaxies, the 
last phase in the metallicity evolution is always associated 
with well after their first pericentre passage and is near or 



1 We thank the referee for pointing out the importance of 
diffusion in this problem. 

2 This is the average distance from the galaxy centre of the 
gas that strongly participates in the interaction-driven inflow. 



after coalescence. This final state of the interaction is true 
whether one considers mergers or simply fly by s. 

3.3. Relationship between Separation and Metallicity 

While the correlation between SFR (at the peak of inten- 
sity) and metallicity dilution is quite remarkable, we expect 
a larger scatter when plotting the dilution as a function of 
the physical distance of galaxies in the pairs, especially for 
those at the smallest separations. This is because the in- 
teracting pair can be in a number of different evolutionary 
states at any given distance, except for perhaps the smallest 
separations. For example, a galaxy at moderate separations 
of several 10s of kpc could be in an initial encounter and be- 
fore first pericentre passage when the circumnuclear metal- 
licity is largely unchanged. Or it could be just after the 
first pericentre passage when its circumnuclear gas could 
be highly diluted (FigJSJ). This causes a range of possible 
metallicity dilutions, from z/zi SO = 1 to z/zi SO ~ 0.5 (FigJSJ 
right panel), with an average value of z/zi SO = 0.75. At the 
smallest separations, say less than 10-20 kpc, but still dis- 
tinguishable systems, galaxies are likely to be either near 
the pericentre passage or near final coalescence. This leads 
to both a strong dilution and, on average, lower circum- 
nuclear metallicity (see FigEJ). At very small separations, 
less than about 4 kpc, galaxies can be near first pericentre 
passage or near coalescence which in either case have rel- 
atively large dilution and relatively low-metallicity, or can 
have coalesced and be strong enriched by supernovae ejecta 
and have relatively high-metallicity. We show this in Fig. [5] 
where in the smallest separation bin, if we include galaxy 
pairs with projected separations of less than 4 kpc we find 
a net average increase in the metallicity (i.e., z >1, with 
an average value of z/zi SO ~ 1.2) while if we exclude such 
systems, we find a net decrease in the circumnuclear metal- 
licity (z/zi SO ~ 0.7). This indicates the predominance of 
the metal enriched mer gers at small projected separations. 

iKewlev et al.l (|2006f ) have shown that galaxy pairs with 
small projected separations (4 kpc h _1 < s < 20 kpc h _1 ) 
and strong bursts of star formation have metallicities lower 
than the comparable field galaxies. Note that in this anal- 
ysis we are considering the physical (3D) distance be- 
tween the two galaxies, and not their projected separation. 
Considering the effects of projection would likely contribute 
to increasing the scatter in the relation between metallicity 
and separation. 

3.4. O and Fe enrichment 

Our simulations do not track the evolution of individual 
elements directly in the code. However, to track their evo- 
lution in an approximate way, we can use some simple ar- 
guments to give an estimate of the amount of iron and a 
elements, like oxygen, released in the ISM by type I and II 
SNe. At a given time, t, during the evolution of a galaxy, 
the rate of oxygen mass released in the ISM by SNe II can 
be estimated by: 

OsNeii(t)[M@/yr] = m 0y SNeiiRsNeiiM*(t) (3) 

where mo y SNeii is the total mass of oxygen ejected aver- 
aged over the initial stellar mass by a type II SNe, RsNeii, 
is the number of SNe II per unit mass of star formation, 
and M*(t) is the mass of stars formed at time t. Similarly, 
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Fig. 3. The duration that the circumnuclear metallicity is below 0.8 (shown in red) and between 0.8 and 1.0 (shown in 
black) relative to an isolated simulated spiral with the same properties as the individual merging galaxies. We chose 0.8 as 
the th r eshold because this is the span of the most extreme metallicity offset observed in nearby mergers (e.g lKewlev et all 
(120061 ): iRupke et al.l (j2008ft ). 



the rate of iron mass ejected on average in the ISM by SNe 
II explosions is given by: 

Fe S Neii(t)[M@/yr] = m Fe , S NeiiRsNeiiM*(t). (4) 

The progenitors of SNe II have very short lifetimes, less 
than a few 10 Myr, so that we can assume that the SNe II 
rate follows closely the star formation rate of the galaxy, 
and so that the release of oxygen and iron by SNe II in the 
ISM is approximately instantaneous 

The stars formed at time t contributes also in enrich- 
ing the ISM with Fe released through type la SNe, which 
are the product of a thermonuclear detonation of a white 
dwarf accreting matter from a companion star. While the 
progenitors of SNe II have very short lifetimes, the time 
delay between the birth of the progenitors of SNe la and 
their explosion can be quite long, of the order of 1 Gyr 
or more d epending on the pro genitors, the binary param- 
eters, etc. (jStroker et al.ll2010D . We will assume delay time 
between the birth of their progenitors and the explosion is 
uniformly distributed between 1 and 4 Gyr after the birth 
of the parent population, thus obtaining that, at any time 
t, the amount of iron released by SNela is: 

FeSNeIa(t)[Mo/yr] = m FeyS NeIaRsNeIaM*(t - t). (5) 

where r is the delay time between the birth of the pro- 
genitors and the detonation of the type la supernova. For 
the supernov ae rates RsNe i a and RsNeii we have used 
the values in IStrolger et al.l ([20101 ) , while the synthesised 
masses m Fe , g AT e Ta 2 rn Fe.s NeTT and m p, s Ne Ta are derived 
from Table 3 in llwamot o et al .1 ([ 19991 ). The amount of oxy- 
gen released by a SNe la is about 0.1 times lower than the 
iron and we will neglect its contribution in the following 
estimations. In Fig]6j we show the rate at which O and 
Fe are ejected into the ISM during some of the simulated 
mergers. We caution the reader that we are showing the 

3 Any possible time delay would result into a time-shift be- 
tween the SFR peak and the oxygen and SNell-Fe peak. 



global values, i.e. evaluated over the whole galaxy. The in- 
tense episodes of star formation are accompanied by in- 
tense peaks of O and Fe release, which coincide with the 
first phases of the encounter, soon after the first pericentre 
passage, and with the merging phase. During the burst of 
star formation, the ISM is enriched mostly in oxygen (the 
OsNeii ejection rate is about a factor 30 greater than that 
of Fes Neil)- If a second burst of star formation occurs in 
the simulation, a concomitant second peak of O and Fe 
release will also take place in the final phases of the inter- 
action causing further relative enrichment of oxygen. 

We note that if the coalescence of the two galaxies hap- 
pens on a short time scale (i.e., less than 1 Gyr after the 
first pericentre passage), there is not enough time to enrich 
the ISM by SNe la. This means that the stellar population 
forming during the final burst forms from gas still rela- 
tively enhanced in a— elements (see FigEJ second panel). 
If, in turn, the merger happens on longer time scales (i.e., 
longer than 1 Gyr), the iron released by SNe la explosions 
has time to lower the [a/Fe] ratio of the ISM, and the 
stellar population formed in the merging phase will be less 
enhanced in a elements (FigEJ fourth panel). 

The O and Fe ejected into the ISM obviously affects the 
ratio between these two elements (FigJTJ). Assuming that 
the initial galaxies have percentages of O and Fe similar to 
the Sun, as star formation proceeds during the interaction, 
the mass ratio of these two elements changes: in particu- 
lar the oxygen abundance increases, until the iron release 
from SNela explosions inverts the trend. Note that the in- 
fluence of the iron release from SNe la on the populations of 
new stars formed during the merger depends on the delay 
time between the birth of the parent population and the 
onset of SNe la - shorter delay times cause a more rapid 
decline in the O/Fe ratio. We caution the reader that in all 
our discussion of the relative abundance of a— elements, we 
are neglecting any possible contribution from SNe la from 
the pre-existing old stellar population. This is of course an 
important simplification of these models which we will in- 
vestigate in future studies. 




Fig. 4. The star formation rate, Mq yr _1 , versus the relative metallicity of the interaction and merger. As in Fig. [TJ 
the metallicity is shown relative to an isolated disk with the same initial conditions as the galaxies used in the merger 
simulations. The blue asterisks indicate times when the star formation is near its peak, while green triangles represent 
galaxies before their peak star formation and red dots represent galaxies that are past their peak of star formation and 
are in the final stages of the merger. In the bottom panels, we show the morphologies of the merger pairs in each region 
of the diagram to indicate the range of morphologies seen in the simulations. 



4. Discussion 

Overall, our results are not surprising. It has been known 
for the last few decades that mergers of gas-rich mas- 
sive galaxi es lead to transformations in the morpholog- 
ical type dBarnes fc Hernauistl 119911 : iBarnesI 1 1 9921 ) and 
drive strong inflows of gas into the circumnuclear region 
(ICasoli et a.1.1 Il99lt iBraine fe Combesl 119931: [Barnes! 120021: 
iBarnes fc Hernauistl 119961 : iDi Matteo et al.l |2007|). These 
inflows are a result of strong tidal torques, dissipation, 
and stro ng asymmetries in the stellar distributi on (e.g., 
bars; see IBarnesI [20021 : [Barnes & Hernquist 1996) leading 
to efficient star formation (|Di Matteo et al.ll2007l ). While 
the infall of gas in mergers has been di scussed within the 
context of simulations and observations ([Casoli et al.lll99ll 



Braine fc Com bes 1993; Barnes 2Q02|; H arnes fc H ernquist 
1996HDL Matteo et al.ll2007l h currently, there has been only 



limited quantitative attempts to address the relationship 
between infall, star formation , and the metallicity (e.g., 
iPerez et al.l 120061 : iRupke et al1l2010D . The trends we find 
in the simulations are fascinating and explain a number of 
important phenomenon that have been observed in inter- 
acting and merging galaxies. Metallicity changes during the 
merger event is a powerful way of testing the timing and 
strength of the gas inflows and ultimately in understanding 
the fate of post-merger/interacting galaxies. 

Overall, we find broad agreement with the results 
of our simulations for the evolution of the metallic- 
ity and what is observed. We find for example that 
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Fig. 5. (Left panel:) The physical 3-dimensional separation of the galaxy pairs and mergers versus their relative circum- 
nuclear metallicity. The black line represents the average relation, for physical distances greater than 4 kpc and binned 
every 10 kpc. The blue line represents the average relation when in addition including mergers with physical separations 
below 4 kpc in the smallest bin (0-10 kpc). (Right panel:) Evolution of two mergers (blue and green curves) and a fly- by 
(red curve) in the metallicity-distance plane. The simulations suggest that even at wide separations there are gas inflows 
which dilute the metallicities of the nuclear regions and that this dilution persists until the pairs are strongly interacting, 
overlapping, and near or after coalescence. During the final merging process, a wide range of metallicities can be observed 
but on average will show a much higher metallicity than if they had not interacted or merged. 




t [Myr] t [Myr] t [Myr] t [Myr] 



Fig. 6. Rate of oxygen and iron released in the ISM by SNe II and SNela explosions (in units of Mq yr _1 ) as a function 
of time two different simulations whose star formation histories are shown in FigJTJ The first and final two panels show 
the ejection rate of oxygen and iron for the star formation history in FigJT] represented by the orange and red curves 
respectively. For the iron release, we separate the contribution from SNela (red curve) and SNe II (grey curve) for each of 
the two simulations. The black arrows in each panel represent the time of the first pericentre passage and final coalescence. 
Since the enrichment is instantaneous (i.e., no delay for the evolutionary times of the type II supernovae, the evolution 
of the star formation and ejecta rate are synchronous. 



the net decrease 
of a dex, which 

in me rging and 

20061: iRupke 



in the metallicity is a few tenths 
is comparable to what is observed 
i nterac ti ng galaxies (iKewlev et al.l 
_.:e et all 120081: iMichel-Dansac et al.l 120081 : 
Alonso. Michel-Dansac. fc Lambas 1 12010ft . In addition, as 
is likely observed, our simulations show that the peak 
of the star formation is close in time to the minimum 



in the circumnuclea r metallicity ([Rupke et al.l [2008; 
Rodri gues et ail 12008ft . The trends could be explained 
entirely by the strength of the gas inflows of relatively 
metal-poor gas from the outer regions of the disks which 
both fuel star formation and dilute the metal content of 
the gas in and around the nucleus. The strength of the 
inflows likely depends both on the initial gas content of the 
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Fig. 7. The evolution of the oxygen over iron mass ratio for two of our galaxy mergers, whose rate of mass release in the 
ISM is shown in FigEl For each panel, three different curves are shown corresponding to three different delay times for 
onset of SNela, r: r =1-4 Gyr (black curve), r = 1-3 Gyr (dark grey) and r = 1-2 Gyr (light grey), as indicated in the 
legends. The black arrows in each panel represent the time of first pericentre passage and coalescence. 
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Fig. 8. Mon t e-carl o simulation of the metallicity distribution predicted by the models compared to the data from 
iRupke et al.l (|2008f ). The up-side down tria ngles are the s imula ted results (see text for details), th e blue stars and red 
circle are galaxy values taken directly from IRupke et al.l ((2008) as described there. Unfortunately, IRupke et al.l (|2008f ) 
did not actually measure the mass of any of the individual galaxies but used a single mass allowing for some scatter (see 
that paper for details). We have approximately reproduced their mass distribution for this comparison. 



disks, how it is distributed, the characteristics of the orbits 
through which the galaxies merge, and the mass ratios 
of the galaxies. Here we concentrated on the modelling 
a range of orbital parameters and a 20% gas fractions. 
Subsequent papers will focus on varying the mass ratios, 
gas fraction, and gas distributions. 

The simulation results suggest that the strongest star 
formation and dilution of the initial nuclear metallicity oc- 
curs around pericentre passage of the merger. This effect 
has been observed in that the star formation at small sep- 
arations exhibit a tail to high Ha equivalent widths (e.g., 



IWoods et all 120061 : iBarton et all l2000l ). At this time, we 
also expect the minimum in the metallicity to occur. This 
expec t relationship has been emphasised in iKewlev et al.l 
(200 6]) based on a stu dy of pairs of galaxies. The sample 
of IKewlev et al.l ([2006) is particularly appropriate for this 
type of study as it is a robust sample of pairs of interact- 
ing galaxies and show generally a modest increase in their 
star formation rates and modest likely dilutions. For exam- 
ple, their most intense starbursts, which show significant 
offsets in metallicity, only have offsets of about 0.2-0.3 dex 
compared to non-interacting galaxies with similar B-band 
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absolute magnitudes. This is the amount we find in our 
simulations. 

We can explain the full range of circumnuclear metal- 
l icities observed in more violent merge rs in the local 
(iRupke et al.l l2008f ) or distant universe ([Rodrigues et al.l 
2008). By violent mergers, we mean in this context galaxies 
which are luminous in the infrared and have star formation 
rates several to almost 100 times that of the Milky Way. 
These are the most appropriate for comparison because, as 
we have shown, as the merger advances, the star formation 
rate will decrease and the metallicity increases reaching val- 
ues that are higher than the initial metallicities. 

To show that our models can reproduce the data, we 
have made a monte-carlo simulation of the metallicity dis- 
tribution we might obtain if we observed a set of galaxies 
near their maximum star formation rate and allow for a 
range in their properties and the intrinsic scatter in the 
mass-met allicity relationship. For such a comparison sam- 
ple, we would naturally draw on the most violent star- 
forming galaxies in the local universe which are major merg- 
ers. This is the point in our simulations, where we get the 
most significant influx of material into the circumnuclear re- 
gions and near the peak of star formation. At this stage, we 
can dilute the circumnuclear gas with relatively metal poor 
gas from the outer disk decreasing the overall metallicity 
by about a factor of 2 and relatively brief time (of-order an 
internal dynamical time of about 10 8 yrs). Thus, in our sim- 
ple monte-carlo approach to the metallicity estimates, we 
allow for a decrease of 0.3±0.1 dex. We make this allowance 
as we only have single metallicity gradient in our all our 
models (—0.1 dex kpc -1 ). However, metallicity gradients 
in spiral galaxies range from about 0.0 to —0.2 dex kpc -1 . 
To allow for this variation, we have assigned a range in the 
metallicity allowed of 0.1 dex. We do not know the nature 
of the progenitors either. They could have a range of initial 
circumnuclear metallicities and the mass metallicity rela- 
tions hip has a significant scatter 0.15 dex ([Tremonti et al.l 
12004ft . Thus we allow for an additional source of scatter of 
0.15 dex to reflect the scatter in the mass- met allicity rela- 
tion. Finally, there is also measuremen t uncertainties which 
are about 0.1 dex ([Rupke et al.ll2008[ ). 

We assume that each source of scatter is Gaussian and 
that each can be linearly combined. This assumption is ap- 
propriate since (at the metallicities observed here) the un- 
certainty in the metallicity is not correlated to the metal- 
licity, our numerical simulations are completely indepen- 
dent of the data, and the scatter in the mass-metallicity 
relationship is intrinsic. One such monte-carlo simulation is 
shown in Fig.[8]but all reproduce the da ta fromlRupke et al.l 
(|2008l ) reasonably well. Unfortunately, IRupke et al.l (|2008l ) 
did not actually estimate masses directly for their sample 
but used a random set of masses generated from a typical 
mass and scatter from the literature where then assigned. 
While we could have used their masses directly which would 
have made the agreement look more spectacular, the real- 
ity is that we have simply reproduced the scatter in the 
metallicity by making a set of simple assumptions. Thus 
this agreement between the data and model should not be 
taken literally but demonstrates that numerical simulations 
of the type run here can reproduce the spread in metallicity 
without making any assumptions beyond suggesting that 
the highest levels of star formation in a merger occurs near 
the peak in the strength of the gas flows something that 



every model over the last decades has shown and can be 
argued on robust physical grounds. 

In addition to investigating the general evolution of the 
metallicity, we also adopted a crude model for the evolu- 
tion of Oxygen and Iron during mergers and interactions. 
We were motivated to make these simplistic estimates by 
the fact that following the evolution of a elements like 
Oxygen can be used to "date" mergers. In analogy with 
studies of the stellar populations of early type galaxies, 
bulges, or in components of the Milky Way such as the 
thick disk and bulge, [O/Fe] or [a/Fe] ratios can be used 
as a chronometer telling us over what timescale the bulk of 
the (luminosity- weighted) stellar popu lation formed (e.g., 
iThomas et aLl '2005: Z occali et a l. 2006). Our hope was that 
even a crude model such as the one we have adopted could 
provide useful information about how observationally we 
might go about dating mergers and interactions. 

How this might work is that if the galaxy is in the early 
phases of merging, we would expect to see both overall 
low- metallicity and depending on the distribution of [a/Fe] 
in initial disks, evidence of what radii the gas came from 
through its abundance ratios. The bulge of our galaxies, at 
least in the stars, is enriched in a elements re l ative to the 
sun dLecureur et al.l 120071 : IZoccali et aD 120081 : iRvde et al.l 
120101 ). Moreover, the disk of our Milky Way appears to 
have an enhancement out to rel atively large radii but this 
ratio depends on the radius (see lAcharova et aT]|2010l and 
references therein). Of course these measurements are for 
the stellar metal abundances, not for the gas phase so the 
comparison is not completely fair. It does however at least 
suggests a possible trend that may be useful in determining 
the radii at which freshly deposited gas came from. At any 
rate, we may therefore expect in the early stages, a mild 
a enhancement in the gas. As the merger proceeds, this 
enhancement will become more pronounced, until the late 
stages where it will begin to lessen as the type la SNe be- 
come more important. An interesting question is how much 
of this enhanced material makes its way into each genera- 
tion of star formation. While in principle this can happen 
quickly as the star formation lasts for several 100 Myrs 
(several crossing times) sufficient for the metals to mix with 
the ambient gas, it is not clear if this actually happens. If 
it does, then we would predict that a majority of the stars 
formed would be a enhanced but may have metallicities 
lower than stars formed from gas before the burst occurred. 
Only after about 1-2 Gyr or so would both the metallic- 
ity increase above what an isolated galaxy would have and 
would the a enhancement begin to lessen. Of course this 
exact time scale is uncertain for both the metallicity and 
a enhance because we do not have a firm understanding of 
the delay time for type la SNe. 

Taken at face value, the results of our simulations 
would suggest that we can make realistic bulges and early- 
type galaxies. It is well known that the bulge of the 
Milky Way is both metal-rich [typic ally about [Fe/Hj~- 
0.3 to 0.3 with a me an of about 0] ( Zoccali et al.l 120 08) 
and i s a enhanced (Lecureur et al.l 2007 : IZoccali et "all 
120081 : IRvde et al.l |2010|). Both of these results are con- 
sistent with our simulations modulo the uncertainty in 
the delay times of type la SNe, and the mass and char- 
acterist ics of pre-ex i sting bulges in the merging galax- 
ies (see iPipino et al.l 12009). This general result appears 
to hold generall y for bulges, namely they are metal-rich 
and a enhanced ([Thomas et al.ll2005l ) and normalizing for 
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their mass difference s, very similar to early type galaxies 
([Thomas et al.ll2005l ). Interestingly, there is some evidence 
for a radial gradient in the average metall icity of the stars 
in the galactic bulge ([Zoccali et al.l l2QQ8h . Depending on 
the initial conditions, this may be consistent with our sim- 
ulations. In the simulations it is likely there is a gradient 
in the star formation with radius as the strongest gas flows 
and most intense star formation is in the deepest part of the 
potential well. However, a complete discussion of the metal- 
licity, a enhancement, and radial gradients in bulges and 
early- type galaxies is beyond the scope of this paper, our 
simulation results are at least consistent with these general 
properties of spheroids. We will address these questions in 
more detail in a subsequent manuscript. 

Finally, these results have important implications for 
our understanding of metallicity evolution of galaxies as 
observed at high-redshift. Since we have found that during 
the most intense phases of star formation, the gas phase 
metallicity can be reduced (diluted) by infalling material 
from the outer regions of the merging galaxies. Mergers and 
interactions are known to play an important r ole in initi- 
ating intense starbursts in the distant universe (lLaw et al 



20071: iNesvadba et~aTl 120081 : lLaw et al.l 12009: E pinat et al 



20091 ). In order to estimate gas phase metallicities of dis- 



tant galaxies, they must hav e high surface brigh tnesses in 
their emission line gas (see iLehnert et al.l I2Q09L and ref- 
erences therein ). It is now reasonably well established 
that the mass-metalli city relation appears to have evolved 
over the last 10 Gyrs (ISavadio et~aLll2QQ5l; lErb et aLll2003 
i Mannucci et aLll2009t iRodrkues et al . 2008; H allidav et ail 
20081: lOuevrel et aLll2009f ). However, combining the possi- 
ble limitations in the observations due to surface bright- 
ness dimming and the importance of mergers in initiat- 
ing intense (high surface brightness) starbursts, it may be 
the determinations of the evolution in metallicity may be 
strongly influenced by dilution. If this is so, we would pre- 
dict that, if the importance of mergers increases with red- 
shift (e.g., Ide Ravel et al.ll2009l ). we should see both a de- 
crease in the average metallicity and an increase in the scat- 
ter in the mass-metallicity relation with increasing redshift. 
While we have no modelled a full range of initial conditions, 
our results here suggest that this increasing scatter would 
be about several tenths of a dex which is less than a factor 
of 2 more than the s catter in the relation for local galaxies 
(|Gallazzi et al.ll2005D . 

5. Conclusions 

In this work, we have studied the dilution and subsequent 
enrichment of the interstellar medium during interactions 
and mergers of disk galaxies of comparable mass. Our mod- 
els, which include star formation and chemical enrichment, 
have made it possible to follow, for the first time, the evo- 
lution of the circumnuclear gas phase metallicity from the 
initial interaction until the coalescence phase, thus signif- 
icantly improving the existing simulations and leading to 
insight between star formation, gas inflows and dilution, 
and metallicity evolution in mergers and interactions. 
Our main results can be summarised as follows: 

— The dilution of the gas initially in circumnuclear regions 
of the merger pairs starts just after the first pericen- 
tre passage, as the disk is destabilised and gas flows 
inwards. For mergers, usually a second dilution peak is 



observed in the final phase of coalescence. A significative 
dilution can take place also in fly- by s, just after their 
close passage, mimicking to some extent what we ob- 
serve in our merger simulations. Thus even fly by s may 
be important for the star formation and evolution of 
galaxies. 

On average, the amplitude of the metallicity dilution 
we find (0.2-0.3 dex ) is in agreement w ith observations 
(jKewlev et al.l 120061 : iRupke et al.ll2008D . 
The strongest correlation we found is between the max- 
imum in the SFR and the strength of the circumnuclear 
dilution - pairs experiencing the strongest bursts also 
show the strongest circumnuclear dilution. This can be 
explained in terms of interaction-driven gas inflows from 
the outer disk regions into the galaxy centre. 
The dilution phase (defined as z/zi SO <0.8) is very 
rapid, half of the sample sustains this high dilution 
phase for less than 5 x 10 8 years. 

The gas inflow causes concomitant enhanced star forma- 
tion; such enhanced star formation releases a significant 
amount of enriched material in the ISM. This enrich- 
ment results in merger remnants which have higher fi- 
nal circumnuclear metallicities than the corresponding 
galaxies evolved in isolation. The same can happen for 
fly- by s which experience enhanced star formation rates. 
Type la and II SNe formed during the interaction release 
O and Fe in the ISM, thus changing the average mass 
ratio of these two elements. In particular, we note that, 
if the merger takes place on short time scales (lower than 
the typical delay time of SNe la, which unfortunately 
is not well constrained), the stellar population formed 
during the last SF burst forms from an a— enriched ISM, 
while, for mergers with longer time scales, SNela can 
have the time to release a sufficient amount of iron in 
the gas. The resulting stellar population has thus lower 
values of [a/Fe] ratios. 

The circumnuclear metallicity, as well as the [a/Fe] 
ratios, can be used as an indicator of the timing of 
merger /interaction state, thus helping in disentangle 
projection effects, i.e., galaxies which appears close in 
projection, but that are not strongly interacting. The 
remnants of these mergers appear qualitatively consis- 
tent with what is known about the bulge of our galaxy, 
other nearby galaxies, and early type galaxies generally. 
Thus, in agreement with many authors, it appears when 
considering chemical evolution that mergers can repro- 
duce early type systems. 

If mergers play a significant role in the galaxies that 
have had their metallicities determined at high-redshift, 
then it is possible that dilution is signficantly affecting 
these estimates. 
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Appendix A: Metallicity evolution in the isolated 
galaxy 

Since the goal of this paper is to investigate the role in- 
teractions and mergers play in the gas phase metallicity 
evolution, it is essential to study the evolution of the gas 
properties of simulated galaxies that did not experience a 
merger or interaction (i.e., an isolated galaxy). This is es- 
sential to distinguish between secular processes that drive 
the evolution of galaxies from the effects that occur dur- 
ing the interaction and merger process. For the isolated 
galaxy simulation, initial metallicity profile (given in Eqn]2]) 
evolves rapidly in the first 500 Myr of evolution, due to the 
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Fig. .2. From top- left to bottom-right: Maps of the gas dis- 
tribution of the model gSb galaxy evolved in isolation. Maps 
are shown from t=0 to t=3 Gyr in steps of 500 Myrs. Each 
map is 40 kpc x 40 kpc in size. 




Fig. .3. From top- left to bottom-right: Maps of the stellar 
distribution of the gSb galaxy model evolved in isolation. 
Maps are shown for every 500 Myr over t=0 to t=3 Gyr. 
Each map is 40 kpc x 40 kpc in size. 

compression of gas into density waves and inflow into the 
central regions (FigOJ). At this time, spiral arms and a bar 
are formed, as clearly shown in the maps of the stellar dis- 
tribution (FigO]). The galaxy models studied in this paper 
were not evolved in isolation before starting the interaction. 
However, as shown in Figs. Al and A2, tidally driven gas 
inflows from the outer disk take place after the pericentre 
passage (with our initial choice of the orbital parameters, 
this corresponds to around 400 Myr from the beginning 
of the simulation) and in the merging phase (between 1 
and 3 Gyr with the exact timing depending on the orbital 
characteristics). By this time - 400 Myrs - the metallic- 
ity profile has already reached an equilibrium configuration 
(see Fig.Al), and the galaxy has relaxed. 



